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We" consider  a  one-dimensional  model  of  the  theimoelastic  response  to 

/N  1  V 

laser  heating  by  a  pulse  of  very  short  duration.  He  treat  the  linear 
equations  of  thermoelasticity  by  perturbation  methods  using  the  pulse  duration 
as  the  perturbation  parameter.  The  perturbation  analysis  separates  the 
problem  into  two  time  regimes,  the  short  time  scale,  which  is  of  the  same 
order  as  the  pulse  duration,  and  the  long  time  scale,  whio*~  is  much  longer 
than  the  pulse  duration.  Convenient  analytical  expressions  are  obtained  for 
the  temperature  change  and  displacement  in  both  time  regimes. 
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SIGNIFICANCE  AND  EXPLANATION 


In  this  pspsr  w  study  the  effects  caused  by  heating  a  material  with  a 
single  short  laser  pulse.  We  consider  only  a  one-dimensional  model  but 
consider  both  temperature  changes  and  elastic  changes,  which  are  coupled 
together.  By  use  of  perturbation  analysis  we  can  separate  the  problem  into 
two  time  regimes,  the  short  time  scale,  which  is  on  the  order  of  the  pulse 
duration,  and  the  long  time  scale,  which  is  much  longer  than  the  pulse 
duration.  Convenient  expressions  are  given  for  evaluating  the  temperature 
changes  and  displacement. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  me,  and  not  with  the  authore  of  this  report. 


THERMOELASTIC  RESPONSE  TO  A  SHORT  LASER  PULSE 
John  C.  Strikwerda  and  Andrew  M.  Scott* 

1.  INTRODUCTION. 

In  this  paper  we  study  the  temperature  change  and  displacements  induced  in  a  solid 
material  by  a  short  laser  pulse.  The  determination  of  the  thermoelastic  response  of  solids 
to  laser  radiation  is  an  important  problem  for  many  applications,  especially  as  laser 
pulses  becosw  increasingly* powerful  and  of  shorter  duration. 

He  apply  perturbation  techniques  to  the  equations  of  linear  thermoelasticity  using  the 
pulse  duration  as  the  perturbation  parameter.  This  approach  has  the  advantage  of  yielding 
relatively  simple  expressions  for  the  temperature  change  and  displacssMint  induced  by  a 
laser  pulse  of  short  duration. 

In  other  work  on  this  problem  Bechtel  (1975)  and  Ready  (1965)  considered  the  heating 
caused  by  laser  pulses  but  did  not  consider  displacsa»nt  effects.  Dunbar  (1981)  considered 
the  effects  of  temperature  on  displacements  but  neglected  the  coupling  effect  of 
displacements  on  temperature.  He  also  included  the  modified  Fourier  law  of  heat 
conduction.  (Dunbar's  results  are  in  error  due  to  an  inappropriate  application  of  his 
initial  conditions.)  All  of  the  above  used  a  Green's  function  approach  to  obtain  explicit, 
though  ccagtlicated,  formulae  for  the  temperature  field.  The  perturbation  method  used  here 
has  the  advantage  of  giving  simpler  expressions  for  the  temperature  and  displacements. 
Moreover,  we  treat  the  full  equations  of  linear  thermoelasticity  in  one-dimension  with  the 
modified  Fourier  law.  The  method  is  limited  by  the  requirement  that  the  pulse  duration  be 
short.  The  method  can  be  extended  to  two  and  three  dimensional  problems  and  to  include 
nonlinear  effects,  however,  the  resulting  equations  will  then  have  to  be  solved  by 
numerical  methods. 
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Nayfeh  (1977)  has  traatad  a  problem  similar  to  tha  ona  considered  in  this  paper  by  the 
method  of  perturbation  expansions.  Ha  considered  an  instantaneous  heat  source  distributed 
on  the  surface  of  the  material.  We  compare  our  results  with  Mayfeh's  in  Section  5. 

The  aquations  describing  the  thermoelastic  response  in  one  dimension  to  laser  heating 

are 


(1.1)  (l  ♦  T  t)K4+  <31  +  2»»>tT0  -  (l  +  tq  ~)a( t/b)e~ °  * 

0  3t  3t  3x3t  tot  3t 

2-  2-  * 

(1.2)  P  -Mr  -  (1  ♦  2u)  “f  +  (31  +  2U)Y  “  -  o 

3t  3x  3x 

where  T  is  the  temperature  change  measured  from  the  Man  state  TQ  and  u  is  the 
displacesMint.  The  solid  material  occupies  the  region  with  x  >  0.  The  normal  stress  is 
given  by 

(1.3)  £  -  (X  +  2y>  —  -  (31  +  2m) YT  . 

3x 

These  aquations  are  those  derived  by  Lord  and  Shulman  (1967)  with  the  addition  of  the  term 
involving  A  which  represents  the  laser  heating.  (He  have  ignored  the  speed  of 
propagation  of  the  laser  pulse  which  is  sufficiently  fast  to  be  regarded  as  infinite.)  The 
parameter  b  is  a  measure  of  the  pulse  width.  The  exponential  damping  of  the  laser 
heating  in  space  is  the  result  of  the  attenuation  of  electromagnetic  radiation  in  a 
conductor.  The  parameter  is  the  relaxation  time  for  the  modified  Fourier  law  for 

which  the  heat  flux  Q  satisfies 

(1.4)  Q  ♦  T.  -  -K  grad  T 

0  3t 

where  K  is  the  thermal  conductivity  of  the  material.  The  other  parameters  in  equations 
(1.1).  (1.2)  and  (1.3)  are  the  Land  elastic  constants  1  and  u.  the  density  p.  the 
specific  heat  at  constant  volume  cy,  and  the  coefficient  of  thermal  expansion  Y>  The 
bars  on  the  several  variables  indicate  that  these  are  dimensional  quantities. 
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Ws  now  non-dinsnsionalizs  ths  equations.  The  natural  length  ecale  for  thin  problea  la 
the  attenuation  length  a  \  so  define  x  «  J/x^  where  is  a-1 .  There  are 

several  poesible  tine  scales  and  we  choose  the  tine  scale  of  heat  conduction  as  the 


reference  tine. 


t  ■  pc  x  -/K  . 
ref  v  ref 


The  perturbation  paraneter  is  to  be  0  which  is  b/tr({,  the  non-dinensional  puls* 
width.  The  relaxation  tine  is  non-dinsnsionalised  as 

T  “  T0/(tref  W  *  V*  * 

The  pulse  nagnitude  is  nornalised  to  have  a  non-dinensional  integral  equal  to  unity. 
Therefore  we  define  ,  by 


/  A(t/b)dt  -  A  -t 


refref 


A(t/b)  -  A^f  -g  A(t/0) 

and  thus 


/  A(s)ds  •  1  . 


The  reference  tenperature  and  diaplaceawnt  are  chosen  as 


Trsf  “  Aref<xref>  ^ 


°r.f  “  *******  *  2‘*>t«f/<‘*r.f>  . 


The  natural  reference  velocity  is 

^r.f  *  Xref/(9tref>  * 

We  thus  obtain  for  the  non- dine nsionallced  equations 


.2-  *2. 


♦  #T  fewft ♦  *  fefc)  -  H  -  o ♦  #T  fe)  h 
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where 


T.  1/2 

fc-)  7  velref  ‘ 


The  non- dimensional  stress  is 


T  m  ± i  1_  *p 

E  "  ta  '  “T  * 
c 


For  many  common  materials  the  non-dimensional  parameters  e  and  d  are  of  order  sero 


in  0,  that  is. 


a  «  c,  d  <<  a 


The  relaxation  parameter  x  plays  an  unimportant  role  in  the  analysis  and  satisfies 

0  <  0x  «  1  . 

The  justification  of  our  treatment  of  these  equations  and  the  non-dt mens Iona ligation  we 
have  used  rests  on  the  above  observations  and  on  the  reasonableness  of  the  results. 

The  boundary  conditions  for  equations  (1.6)  and  (1.7)  are  the  conditions  of  no 


radiative  heat  loss 

(1.10) 

and  no  stress  at  the  wall 

(1.11) 


£  “  0  «  x  -  0  , 


8u  1  _  _  _ 

m-*2*'o  at  *  - 0 


The  initial  conditions  are  that  T(t,x)  and  u(t,x)  and  all  their  derivatives  vanish  as 
t  approaches  negative  infinity. 

He  will  solve  equations  (1.6)-(1.11)  using  perturbation  expansions  with  0  as  the 
perturbation  parameter.  Since  0  is  on  the  order  of  10-10  for  many  problems,  excellent 
results  can  be  obtained  with  only  one  term. 

The  paper  is  organised  in  the  following  manner.  In  section  2  we  derive  the 
perturbation  expansion  for  the  short  time  scale,  that  is  times  on  the  order  of  the  pulse 
width,  and  in  section  3  we  consider  the  long  time  scale  behavior.  In  section  4  we  solve 
for  the  boundary  layer  correction  at  the  material  surface  which  is  needed  to  complete  the 
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2.  THE  SOLUTION  ON  THE  SHORT  TINE  SOLE 


2  *S  .  ",  „ 

la.Jl 

3t? 

^r*»r-0 

with 

3Tt 

*“1  -2 

(2.4) 

•jj-  *  0  and 

C  1,-0  at  x-0 

and 

atk  '  at1* 

(2.5) 

0  as  t  ♦  for  k  >  0  . 

We  Htk  solutions  to  those  equations  *n  the  fora  of  power  series  in  ths  perturbation 
parameter  0  i.e. 

T  -  T®<t1#x)  ♦  0r](t1#x)  +  02T*(t1,x)  +  ... 

u  -  Ujttj.x)  ♦  0u’(t1#x)  +  02uj(t1,x)  ■*■  ...  . 

Substituting  these  expansions  into  equations  (2.2)  and  (2.3)  and  equating  like  powers  of 
0  gives  the  following  system  of  equations 


(2.6) 

(2.7) 


(T  if: +  -  A<V*"X ♦  d2  - 0 


»2u*  2  32U^  ^ 

3t2  3x2  +  ** 


-  0 ,  for  k  >  0 


and 
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equation  (2.6)  can  ba  intagrated  to 


3T 


-  , 


./T 


for  naa  function  C(x)  Tha  initial  conditions  (2.6)  iaply  that  C  ia  idantically  saro. 
Another  integration  in  t  and  using  (2.6)  again  gives 


(2.9) 


1  Ju 

°(t1#x)  -  a"*  /  A(s)ds  -  d2  -gp  (t,,x> 


Diffarantiating  (2.9)  with  raapact  to  x  and  substituting  in  (2.7)  for  k  ■  0  gives 

t. 


(2.10) 


with 


.2  0 
a  u 


,2  0 


— r1  -  (e2  ♦  d2)  — r1  -  •“*  /  A(s)ds 


at: 


ax' 


(2.11) 


ax 


(e2  ♦  d2)  -• 


/  A(s)ds  at  x  -  0  . 


The  aolution  to  (2.10)  and  (2.11)  satisfying  the  initial  condition  (2.5)  is 

t. 


(2.12) 


(t^.x)  •  e”x  /  a”2 (cosh  a(t^  -  s)  -  1)  A(s)ds 


r  -2 

-  /  a  sinh  a(tf  -  s)  A(s  -  x/a)ds 


a  ■  (o2  ♦  d2)1^2.  rron  (2.9)  wa  obtain 


*1 

(2.13)  T°(t1#x)  ■  c2/a2e_x  /  A(s)ds  ♦  2/a2  /  co*v  j(t1-e)(a-xA(s)  -  A(s-x/a))ds 
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-8' 


f (r)  -  /  *“MA(s)ds  , 

-•  r 

-  r.-‘|r-*|A(.)d. . 

Not*  that  the  twptrttur*  change  associated  with  the  compression  wav*  ia  first  positive 
than  negative  as  the  wav*  passes  any  point. 

The  functions  T*  do  not  satisfy  th*  boundary  condition  (1.10).  A  boundary-layer 
expansion  is  n**d*d  to  satisfy  this  boundary  condition  and  this  expansion  is  the  topic  of 


section  4, 


In  this  section  we  consider  the  theraoelastic  response  on  the  long  time  scale,  that  is 
for  times  much  larger  than  the  laser  pulse  width.  The  equations  to  be  solved  are  (1.6)  and 
(1.7)  which  we  rewrite  as 


(3.1)  If  -  +  d2  |jrr  +  8t  +  d2  -^-z)  -  \  (A(t/6)  +  t  A’(t/W)e 

3t  3x2  3x3t  3t2  3x3t2  6 


r  32T 


33u 


(3.2) 


2 

-  c 


0  , 


with  the  initial  conditions  that  T  and  u  vanish  as  t  approaches  negative  infinity, 
and  the  boundary  conditions  (1.10)  and  (1.11).  We  will  obtain  the  solution  as  a 
perturbation  expansion  in  the  variable  6.  However,  since  the  right-hand  side  of  (3.1)  has 
a  singularity  at  6  -  0,  the  terms  of  the  perturbation  expansion  will  be  distributions  and 
we  must  consider  the  weak  form  of  (3.1)  with  regard  to  derivatives  in  t. 

First,  let  us  consider  equation  (3.2).  Writing 


and 


T(t,x)  ~  l  0M(t,x) 
j-0 


we  obtain 
(3.3) 


u(t,x) 


-  I  PJuJ(t,x) 
3-0 


2  3V 

c 


♦ 


2  j-2 

•I 


3  u 


3  -  0,1 


—2  —1 

where  u  «  u  ■  0. 


When  j  is  0  and  1,  we  have,  by  integrating, 


which  satisfy  both  tha  no  atraaa  boundary  condition  (1.11)  at  x  •  0,  and  tha  dacay  of 
u  and  T  as  x  ♦  •.  Thus  to  within  0( B2)  thara  is  no  straws  on  tha  long  tins  seals. 
Tha  weak  for*  of  (3.1)  which  wa  consider  ia 


m  2 

/  (-r(t)T  -  4(t)  -2-|  -  d2#'(t)  *  0Tf  (t)(T  «•  d2  ~))«lt 

••  3x 


-4  /  (*(t)  -  0T*'(t))A(t/8)dt  a" 


for  all  +  e  c"(R).  To  continue  wa  auat  expand  tha  right  hand  aids  of  (3.5)  aa  a  aariaa  in 
0.  Making  tha  change  of  variable  t  -  Ba  wa  have 


4  /  (*(t)  -  0T*'(t)>A(t/0)dt 

P  — m 

m 

-  /  (♦(Be)  -  0T«>(0a))A(s)ds  . 


Mow  for  any  positive  integer  M  wa  will  obtain  an  expansion  of  tha  above  integral  with 
error  0(  0**).  First  choose  R  so  that 

I  /  skA(s)ds|  <  rf*"k  for  k  -  0.....M  . 

|a|>R 

Than,  by  tha  Taylor  series  expansion  of  +(0s). 


•  R 

/  +(0s)A(s)ds  -  /  +(0s)A(s)ds  4-  /  ♦( 6s)A(s)ds 

•  -R  |s|>R 

M-1  R  ,  m  , 

-  I  I  ♦'3,(0)(Bs):JA(s)ds 
j-0  -R  3 

R  i 

♦  P  I  /  ♦<”) ( Bs' )(s  -  s' J^’ds'ds  +  /  «0s)A(s)ds 

-R  1,1  o  |s|>r 


M-1  A(j)|n.  •  ,  „ 

-  I  B3  *  .  <0>  /  83A(s)ds  «■  0<  ^*)  «♦■„ 

j-0  31  —  '* 
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where  l$l  *  sup  sup  '(s)|.  For  convenience  define 

#M  s  0<  j<M 

00 

(3.6)  “j  “  jT  ^  sjA(a)da  j-0,1,... 

with  a_«  -  0.  By  the  normalisation  of  A(t)  we  have  that  aQ  ia  1. 


Therefore, 


and  similarly 


00 

4  /  ♦(t)x(|)at  -  l  eV 4(3> co)  ♦  o< «?*)»♦».  t 

p  -•  p  j-0  3  ' 


00 

/  4'(t)A(t/B)dt  -  J  bV  ,VJ,<0>  +  0(fifl)IM  „ 

—  j-1  3-1  »" 


Thus  the  right-hand  aide  of  (3.5)  ia 


M~1 

I  B3<a.  -  Ta.  ,)4<j>(0)  ♦  0(«f*)l*l  „  . 
j-0  3  3  1 

Now  expanding  the  ^eft-hand  aide  of  (3.5)  in  powers  of  B  we  have  for  j-0 


-  .  .2.0  „  .  0 

/  -4*t°  -  4  -  4  4’  ■§-  dt  -  4(0 )« 

~m  ax 


and  using  (3*4)  we  obtain 


/  ”4'(l  ♦  %)t°  -  4  dt  -  4(0)e"x 


Since  the  right-hand  aide  of  this  equation  depends  on  4(t)  only  at  t  -  0  we  have 
that  ia  differentiable  for  all  t  other  than  t  -  0.  Integrating  by  parts  then  givea 
us 

/  4<t){(l  +  %)  |f-  -  ^-?-}dt  +  4(0)[T°(0+,x)  -  T°(0-,x))  -  4(0)e”x  . 
c  3x* 

Thus  T°  is  a  solution  to  the  heat  equation  for  t  >  0  and  t  <  0,  with  the  homogeneous 
boundary  condition  (1.10).  Since  T°  vanishes  aa  t  ♦  —  the  uniqueness  of  the  initial 
boundary  value  problem  implies  that  T°  is  identically  zero  for  t  <  0.  Therefore  for 
t  >  0,  T°  satisfies 


•». 
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with  T  (0,x)  -  k*  and  jt-  (t,0)  -  0  where  k  -  cV(c*  ♦  d2). 

Mot*  that  tha  initial  condition  for  T°  ia  tha  same  aa  tha  limit  of  aa  t, 
baoomas  infinite  for  tha  abort  aoala  (2.15).  similarly  for  u®,  from  (3.4)  and  tha 
initial  value  for  T°  wa  have 

(3.8)  u°(0,x)  »  -  -r^ — r  a"*  > 

c  ad* 

which  agree*  with  (2.14).  Hi*  explicit  eolation  to  (3.7)  ia 
(3.8)  T°  (t,x)  -  \  (.-<*-y)2/4kt  +  €-(xay)2/4ktu-yay 

-  |  *kt(*“*(1  -  erf(/kt  -  x//tkt)>  ♦  a*(1  -  arf (^t  a  x//Sct))>  , 

and  ao 

a  1  «  • 

(3.10)  u  (t.x)  -  ■  (-1  a  ^  /  (arf((x  -  y)//4kt)  a  arf((x  a  y)//3kt) )*“ydy )  . 

o*  a  d*  0 

From  (3.4)  it  ia  aaan  that  u°  la  detarminad  only  up  to  a  function  of  t.  Xn 
deriving  (3.10)  we  have  uaad  the  condition  that  u®  vanlahaa  aa  x  baoomaa  larger.  From 
(3.10)  wa  aa*  that 

u°(t,x)  a  -i/(c2  a  d2)  aa  t  *  •  . 

This  non-sero  limit  for  the  displacement  la  tha  result  of  naglaoting  any  raatraining  forces 
on  the  material  which  would  restore  it  to  its  original  position.  Tha  aquation  for  T®, 
using  (3.4),  is 

/  (-4'T®  -  4k  a  4«TT°)dt  -  4* (O)k(a  -  •»„)*"*  . 

-•  3x2  10 

As  before  T®  is  saro  for  t  <  0.  T®(t,x)  can  be  written  as  tha  sum  of  a  distribution 

and  a  smooth  function  T®(t,x).  Ns  have 


T®(t,x)  -  -ka,«(t)*"x  a  T*(t,x) 


where  6(t)  is  ths  usual  Dirac  delta  function  and  T  (t,x)  satisfies 


tjV  .  o 

H  »x2  at2 


with  T1  (0,x)  “  -k(a^k  ♦  t)s  x  and  -jj-  (t,0)  -  0.  Ths  aquation  for  u^(t,x)  from 

is 

m 

u’tt/X)  -  — r  ka1«<t>s-x - j  /  T^t.yjdy  . 

c  *  ex 

In  a  similar  way  ons  may  derive  expressions  for  T2  and  u2,  however,  for  most 


.4) 


applications  the  first  two  terms  should  suffice 


* 


4.  TIB  BOUNDARY  COMSCTXON  Oil  THB  SHORT  TUB  SC Aid. 

Xn  this  section  we  calculate  the  correction  teres  which  are  needed  to  satisfy  the 
radiation  boundary  condition  (1.10)  for  the  short  tins  scale.  As  in  section  2  the  ties 
variable  is  tf>  l.e.  t/6,  and  we  write  the  teasers ture  function  as 

(4.D  Tit^n.e)  -  T^t^x.o)  ♦  i1/2  T1(t1.x*r1/2,a)  . 

Setting  x^  ■  xfl”1^2,  we  have  by  (1.10) 


(4.2) 


**1  *1 

(tj.0,0)  ♦  (tj  ,0, 0)  ■  0  . 


satisfies  the  sane  Initial  conditions  as  Tt,  that  is  t(  and  all  its 
derivatives  vanish  as  tt  tends  to  negative  infinity. 


To  obtain  neaningful  equations  for  ^  and  u^  we  set 


(4.3) 


u(t1  fX,  0)  -  u1  (t|  fXf  0)  ♦  0u1(t1,x1,0) 


and,  free  (2.2)  and  (2.3),  the  equations  for  and  are 


(4.4) 


2*  *  2* 

•V  «>  >  *■< 

rqs; 


*2*. 


(4.S) 


The  aystee  (4.4)  and  (4.3)  with  the  boundary  condition  (4.2)  and  the  initial 


conditions  determine  Tj  and  u^ •  He  again  solve  these  equations  using  perturbation 
expansions  in  0.  He  write 


*.  (t,  fx, ,  0)  -  l  (^(t.fX.) 
fc-0 

(4.6) 

u1(t.,x1,0)  -  l  ^u^(t1#xt)  . 

k-0  1,1 
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•re 


*0  “0 

The  equation*  for  T  end  u 


and 


,2*0  *0  ,  ,2*0  ,  ,3  *0  ,2*0 

9  T  9T  2  9  u  2  9  u  9  1 

T  — T~  *  TT  *  d  ZT  Zr  *  Td  — H - T 

9t2  9t1  9t19x1  9t29x  *** 


**  'sr0 ' 


Integrating  this  laet  equation,  we  have 

*0 

,  9u.  «„ 

2  10 

_c  _-+Ii.0  , 

which  ehows  that  there  ie  no  correction  to  the  atreaa  to  the  firat  order  in  6.  We  then 

have  2  ,q  ,q  2 ,  q 

3  T.  9T.  3  T, 

(4.7)  T_TL  +  _-L.  H— -1-  0 

9tJ  *1  9xJ 

2  2  2  *0 

where,  aa  in  section  3,  k  -  c  /(c  ♦  d  ).  The  boundary  condition  for  T1  is,  from  (4.2) 
and  (2.13), 


(4.S) 


3*.  1  -at.  1 

■s-1  (t  ,0)  -  k  /  A(a)  ♦  (1  -  k)e  /  ea*A(s)da  -  B  (t. 

.«  1 


The  solution  to  (4*7)  and  (4*8)  is 


.WC1  -(t  -s)/2T  ,  x,  2  1/2  „ 

'?<VV  -  -c ,  /  •  I0(77  ((t,  -  a)2  -  (~)  )  )  B°(a)da  . 

-•  1 


where  c,  ie  /  k/x  and  IQ  is  the  Modified  Bessel  function  of  order  sero.  The  foraula 
for  the  displacement  correction  la 

c.  2  tl"x1/c1  -(t  -s)/2i 

u1(t1,*1)  ■  (— )  /  e  0<t1  -  s,x^/c)B°(s)ds 

where 
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We  now  preaent  mm  explicit  results  of  ths  preceding  analysis.  We  consider  two  pulse 
shapes i  the  first  pulse  has  the  shape  of  a  Gaussian 


(3.1)  A(s)  -  s"VV* 

and  the  second  has  an  exponential  decay 


(5.2) 


A(s)  - 


0 


-a 


s  <  0  , 
s  >  0  . 


The  analysis  of  sections  2  and  3  shows  that  one  may  take  A(s)  to  be  any  function  with 
unit  integral  satisfying 

1)  /  skA(s)ds  <  •  k  -  0,1,... 


2)  e4t  /  e“**A(s)ds  *0  as  t  ♦  — , 
where  a  -  (c2  ♦  d2)1/2. 

Graphs  of  T®  and  u®,  the  short  tins  scale  temperature  change  and  displacement,  are 
shown  in  figures  1  to  5  for  several  values  of  x  as  functions  of  t, .  The  values  of 
T®  and  nj  are  computed  using  equations  (2.12)  and  (2.13)  for  x  -  0,1,2,  etc.  The 
values  of  t^  run  up  to  10.0  in  each  figure. 

figures  la  and  1b  show  the  temperature  change  and  displacaswnt  for  the  Gaussian  pulse 

(3.1)  for  values  of  c  ■  1.0  and  d  “  1.0.  Figures  2a  and  2b  show  similar  results  for  the 
exponential  pulse  (3.2)  with  c  •  1.0  and  d  «  1.0.  Comparing  figures  1  and  2  one  can  see 
the  effect  on  the  temperature  and  displacement  of  the  different  rates  of  energy  depoeition 
due  to  the  different  pulse  shapes.  The  exponential  pulse  (5.2)  deposits  energy  more  slowly 
than  does  the  Gaussian  pulse  (5.1)  and  therefore  the  temperature  and  displacement  approach 
their  limiting  values  more  slowly.  Figures  3a  and  3b  show  results  for  the  Gaussian  pulse 

(5.1)  for  c  “  0.7  and  d  “  1.0.  Figures  4a  and  4b  display  the  temperature  change  and 


iJ.'JJJWl 


.ms  than  10,  c-0.7  and  d»1.0 


I 


3 


h  ■** 


C-j 


displacement  for  the  Gaussian  pulse  (5.1)  with  c  «  0.7  and  d  *  0.7.  Figure  5  shows  the 
temperature  change  for  the  Gaussian  pulse  (5.1)  with  c  »  1.0  and  d  ”  0.7.  The 


displacement  corresponding  to  figure  5  is  the  same  as  figure  3b.  In  each  case  the  limit  of 


22  2  -1  -x 

the  temperature  for  large  t,  is  c  (c  +  d  )  e  and  the  limit  of  the  displacement  is 


/  2  .2-1  -x 

—  (c  ♦  d  )  e  . 


One  clearly  sees  in  these  figures  the  beginning  of  the  compression  wave  given  by 
(2.16)  and  (2.17).  The  negative  temperatures  associated  with  the  compression  wave  are  a 
consequence  of  the  strong  thermoelastic  coupling  and  comparative  insignificance  of 
diffusion  on  the  shoii  time  scale.  Similar  behavior  is  not  seen  in  Nayfeh's  results 
(Nayfeh  1977)  since  in  that  work  the  thermoelastic  coupling  coefficient,  i.e.  d,  is  used 
as  the  perturbation  parameter. 

The  different  figures  also  show  the  effect  of  the  parameters  c  and  d  on  the 
relative  magnitudes  of  temperature  and  displacement  of  the  compression  wave  as  compared  to 
the  temperature  and  displacement  on  the  surface.  Further,  figures  1,  2,  and  4  illustrate 
the  decrease  in  speed  of  the  compression  wave  as  c  and  d  are  decreased.  The 
temperature  on  the  surface  is  independent  of  c  and  d  except  through  a  multiplicative 
factor,  therefore  because  of  our  scaling,  the  temperature  curves  for  x  «  0  are  identical 
in  figures  1,  2,  4  and  5. 

Figure  6  shows  the  long  term  behavior  of  the  temperature  for  values  of  c  and  d  of 
1.0,  hence  k  is  0.5.  Equation  (3.9)  is  plotted  for  values  of  x  of  0,  1.0,  2.0,  etc. 
and  t  ranges  up  to  40. 

The  results  displayed  in  these  figures  agree  qualitatively  with  the  results  of  Bechtel 
(1975)  and  Ready  (1965).  Moreover,  the  final  results  are,  we  believe,  easier  to  compute 
and  more  informative  than  previous  results  since  the  perturbation  analysis  emphasises  the 
dominant  effects. 

The  original  system  of  equations  (1.1)  and  (1.2)  has  two  propagation  speeds  associated 
with  it  when  the  relaxation  time  is  non-sero.  In  our  analysis  these  show  up  on  two 


different  scales.  On  short  time  scale,  as  discussed  in  section  2,  the  effects  propagate 


2  2  1/2 

with  the  non-dimensional  velocity  a  ■  (c  +  d  )  ,  and  in  the  boundary  correction  for 


-24- 


L>'-. 


the  short  time  discussed  in  section  4  the  thermal  disturbancee  propagate  with  speed 
—1/2  2  1/2 

c1 8  -  c(a  T0)  .  Moreover,  the  disturbances  propagating  with  the  faster  speed 

-1/2 

Cj0  die  out  rapidly  and  those  propagating  with  the  slower  speed  do  not  die  out. 

These  results  agree  with  those  of  Nayfeh  (1977)  and  are  compatible  with  the  results  of 
Norwood  and  Warren  (1969).  on  the  long  time  scale  of  section  3,  there  is  no  propagation 
such  since  diffusion  is  the  doaii nant  phenomena. 
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